arXiv:1508.05426v2 [astro-ph.HE] 2 Dec 2015 


Accepted for Publication in The Astrophysical Journal 

Preprint typeset using DTjtX style emulateapj v. 5/2/11 


FERMI- LAT OBSERVATIONS OF SUPERNOVA REMNANT G5.7-0.1, BELIEVED TO BE 
INTERACTING WITH MOLECULAR CLOUDS 

Timothy Joubert 1 - 4 , Daniel Castro 2 , Patrick Slane 1 , Joseph Gelfand 3 
Accepted for Publication in The Astrophysical Journal 

ABSTRACT 

This work reports on the detection of 7-ray emission coincident with the supernova remnant (SNR) G5.7- 
0.1 using data collected by the Large Area Telescope aboard the Fermi Gamma-ray Space Telescope. The 
SNR is believed to be interacting with molecular clouds, based on 1720 MHz hydroxyl (OH) maser emission 
observations in its direction. This interaction is expected to provide targets for the production of 7-ray emission 
from 7r°-decay. A 7-ray source was observed in the direction of SNR G5.7-0.1, positioned nearby the bright 7- 
ray source SNR W28. We model the emission from radio to 7-ray energies using a one-zone model. Following 
consideration of both 7r°-decay and leptonically dominated emission scenarios for the MeV-TeV source, we 
conclude that a considerable component of the 7-ray emission must originate from the 7r°-decay channel. 
Finally, constraints were placed on the reported ambiguity of the SNR distance through X-ray column density 
measurements made using XMM-Newton observations. We conclude G5.7-0.1 is a significant 7-ray source 
positioned at a distance of ~ 3 kpc with luminosity in the 0.1-100 GeV range of L 7 ss 7.4 x 10 34 erg/s. 

Subject headings: acceleration of particles — cosmic rays — gamma rays: ISM — ISM: individual ( G5.7-0.1) 
— ISM: supernova remnants 


1. INTRODUCTION 

Observational evidence supporting the diffusive shock ac¬ 
celeration (DSA) of cosmic rays (CR) has steadily increased 
over the last few decades. Current evidence includes obser¬ 
vations of nonthermal X-ray emission from shell-type SNRs, 
believed to be the result of synchrotron radiation from shock 
accelerated electrons fLazendic et al.[2004|[Vink|2012j ). Also, 
very high energy (VHE) 7-ray emission has been observed 
in the direction of a significant number of SNRs, including 
SN 1006, RX J0852.0-4 622, and IC 443 using ground based 
Cherenkov telescopes (fAharonian et al. 2007 Acciari et al. 
|2009| |Berezhko et al.)2009| ). Both the X-ray and VHE 7-ray 
cases support a nonthermal shock acceleration mechanism, 
one that would result in acceleration of both hadrons and elec¬ 
trons, which can be termed hadronic and leptonic accelera¬ 
tion respectively. There are several emission mechanisms that 
result from the acceleration of electrons including nonther¬ 
mal bremsstrahlung emission and inverse Compton scattering 
(IC) of photons by accelerated electrons. Hadronic accelera¬ 
tion allows for a separate 7-ray emission mechanism caused 
by interactions between the accelerated hadrons and ambient 
material ( |Ellison et al.|2010] l. 

The hadronic mechanism should theoretically produce 7- 
ray emission through the decay of neutral pions, 7r° —> 77, 
produced by proton-proton, generally nuclear-nuclear, inter- 
actions of the shock acce lerated ejecta and ambient material 
( Ackermann et al.]|2013| l. This material often exists in the 
form of molecular clouds (MC) adjacent to the SNR ([Drury 
et al .|l994 Claussen et al. |l997j ). Acceleration of cosmic-ray 
protons into these clouds, or re-acceleration of ambient cos- 
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mic rays, is believed to enhance the 7r°-decay 7-ray emission 
signature due to more frequent pion producing interactions 


because of the higher proton density relative to the ISM (Aha- 

ronian et al. 19941 Drury et al. 1994 

|Uchiyama et al. 2010 

Tang & Chevalier|2014| Lee et al. 20 

15). Therefore, observa- 


SNR-MC systems, provide a particularly interesting astro- 
physical laboratory for studying the signatures of a hadronic 
7-ray production mechanism (Wardle & Yusef-Zadeh|2002l. 


The theorized peak of the resulting photon spectrum 
is well positioned within the detectable energy range of 
the Fermi Gamma-ray Space Telescope Large Area Telescope 
(Fermi'-LAT), and is characterized by a sharp rise between 


- 70—200 MeV (Kamae et al.[2006; Reynolds[2008). Fermi- 
LAT has made significant improvements in sensitivity and res¬ 
olution over previous instruments covering this energy range, 
most notably the Energetic Gamma-Ray Experiment Tele¬ 
scope (EGRET). This makes Fermi-LXT study of SNR-MC 
systems very promising. Further, recent Fermi- LAT observa¬ 
tions of SNR-MC systems have corroborated theory, estab¬ 
lishing these interactions as significant 7-ray sources, includ¬ 
ing correlation studies performed by Hewitt et al. (20091. In 


fact, SNR-MC systems are among the most luminous SNRs 


in 7-rays (jAbdo et al.|2009c Thompson et al. 20111. 
While there are contrasting theories regarding the 


garding the specific 


acceleration mechanism, evidence supporting a hadronic ori¬ 
gin of observed 7-ray emission has recently strengthened. 
Observation of the SNR-MC systems IC443 and W44 by |Ack-| 


ermann et al. 


(2013) clearly defined the 7r°-decay signature 
down to 60 MeVTeffectively demonstrating the hadronically 
dominated origin of the observed 7-ray emission in the direc¬ 
tion of each SNR. Additional support for the connection of 
SNR-MC systems to hadronic acceleration comes from ob¬ 
servations of the SNRs W41, MSH 17-39 and G337.0-0.1 


by Castro et al. (2013). This work compares three unique 
SNR-MC systems and the closeness of fit of two sophisti¬ 
cated 7-ray emission models to the broadband emission spec¬ 
trum of each. The study reports that a model assuming a pre- 
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dominantly hadronic emission mechanism was in close agree¬ 
ment for the spectra of each SNR. Similar support is pro¬ 
vided by the Fermi -LAT study of Kes 17 by Gelfand et al. 
( |2013| ) wherein a leptonically dominated emission scenario 
for the observed 7 -ray emission was ruled out due to the high 
required electron-proton ratio as compared to local observa- 
tions of the ratio. Th e Fermi -LAT study of SNR Kes 79 by 
Auchettl et al. [ (12014b also supported a hadronic origin to the 
7 -ray emission based on both a high required electron-proton 
ratio compared to local observations, and also that the total 
electron energy required for a leptonically dominated emis¬ 
sion origin would exceed the canonical SN explosion energy. 

Evidence of SNR-MC interactions have been observed 
through focus on various emission lines including the detec¬ 
tion of O H(1720 MHz) mas er emission, which was first ob¬ 
served by |Goss & Robinsonj ( |l968) l toward the SNRs W28 and 
W44. This emission has become a defining indicator of SNR- 
MC interaction. While other emission lines may indicate the 
existence of an MC within the field-of-view (FoV) of an SNR 
observation, the precise position of the MC is not easily de¬ 
termined. The OH(1720 MHz) maser emission line however 
is only emitted by dense gas (~ 10 5 cm' 3 ) following the pas¬ 
sage of non-dissociative shocks (|Lockett et al.||1999). This 
makes OH maser observations adjacent or interior to SNR ra¬ 
dio and X-ray contours indicative of shock induced accelera¬ 
tion of cosmic ray hadrons, and potential sites of gamma-ray 
emission from these hadrons colliding with ambient material 
(Claussen et al. 1997*|Koralesky et al.|l998) . 

The SNR W28 has been observed extensively at energies 
ranging from the radio to the very high energy (VHE) 7 -ray 
regime. Within the field-of-view of W28 is the subject of this 
paper, SNR G5.7-0.1 for which very little direct study has 
been conducted at any wavelength. G5.7-0.1 was discovered 
by Brogan et al.] ( 2006) using observations of the the region 
taken with the VLA. Due to the SNR’s relatively low radio 
brightness, G5.7-0.1 was not definitively considered an SNR 
with a radio shell, and was not initially included in Green’s 
catalog (Gree np009 1. The 2006 detection briefly established 
the existence of an SNR candidate named G5.71-0.08 which 
was then changed to G5.7-0.0 following additional studies of 
the region by |Hewitt et al. I ( |2009] >, using VLA and GBT ob¬ 
servations. Data from the Multi-Array Galactic Plane Imag¬ 
ing Survey (MAGPIS) ( Helfand et al.|2006) was also used to 
study the spatial characteristics of it’s partial radio shell. This 
study confirmed the SNR categorization of G5.7-0.1. 

The SNR is believed to be interacting with neighboring 
molecular clouds, as evidenced by OH(1720 MHz) maser 
emission observed coincident with the radio center of the 
G5.7-0.1 ( |Hewitt & Yusef-Zadehl2009] l and also several stud- 
ies of cloud interaction in the region of SNR W28 (Claussen 


et al. 1997 Arikawa et al.jl999 Velazquez et al. 7) 

low radio brightness of G5.7-0.1 has made previou 


02). The 

brightness of G5.7-0.1 has made previous attempts 

at setting constraints on the distance to the SNR difficult. Ve¬ 
locity measurements taken from OH(1720 MHz) maser obser¬ 
vation set a kinematic velocity of 12.8 km s _l . Due to veloc- 
ity ambiguity which could not be resolved, Hewitt & Yusef- 


Zadeh (2009) successfully set kinematic distance constraints 
at either 3.1 or 13.7 kpc. The maser observation establishes 
G5.7-0.1 as a strong candidate site for shock induced cosmic- 
ray acceleration. In addition to maser detection, Aharonian 
|et al.| (j2008|) detected the VHE 7 -ray source HES 
240C inside the radio extent of the SNR while studying the 
W28 region, and also claimed a hadronic origin. 


This paper will study the 7 -ray emission in the direction 
of G5.7-0.1 and seek to establish constraints on the nature of 
the production mechanisms by modeling the broadband spec¬ 
trum. We will also use comparisons of different best-fit model 
parameters and XMM-Newton X-ray analysis of the region to 
set constraints on the distance to the SNR. A description of 
how the Fermi-LAI data were analyzed and the corresponding 
observation will be given in Section [2] followed by a descrip¬ 
tion of these results and their modeling in Section [3] Finally 
a summary of the results is given in Section[4] 


2. OBSERVATIONS AND DATA ANALYSIS 

The Fermi -LAT is an e~e + pair conversion telescope with 
a field-of-view (FoV) of ~ 2.3 sr, and an energy dependent 
point-spread function (PSF) that improves at higher energy. 
The Fermi-LAI is capable of detecting photons with energies 
from 20 MeV to 300 GeV. Further details regarding the in¬ 
strument and its calibration can be found in lAtwood et al.l 
(2009). This work analyzes Fermi-LAI data spanning nearly 
63 months of measurement, from 4 August, 2008 until 27 
September, 2013. Only events with the Pass 7 V 6 Source 
class designation were used during this analysis, which sig¬ 
nificantly reduces the residual background rate (Ackermann 
|et al.|[ 20 T 2 l and references therein). In addition to the event 
class designation, analyzed events had reconstructed zenith 
angle values smaller than 100 ° to reduce the inclusion of ter¬ 
restrial 7 -rays ( Abdo et al.|2009b) . All analyzed events came 
from a circular region, centerecTat the radio center of G5.7- 
0.1, with a radius of 25°. The instrument response functions 
(IRFs) used in the analysis correspond to the "Pass7 version 
6 " versions developed using in-flight data (Rando 2009; Ack- 
|ermann et al. 12012 ) . 

Data analysis was performed using the Fermi Science Tools 
v9r31pl |^| Maximum likelihood analysis was conducted on 
the data set using the gtlike tool, (Mattox et al.|1996(l, to study 
the morphological and spectral properties of the SNR. The 
emission models used by the fitting routines include a Galac¬ 
tic diffuse component describing the interactions of cosmic 
rays with the ISM and cosmic-ray electron interaction with 
interstellar radiation, and an isotropic component describing 
the instrumental and extragalactic diffuse backgrounds. The 
mapcube modeling the Milky Way 7 -ray emission was the 
file gal_2yearp7v6_trim_v0 .fits, with the isotropic 
emission modeled by the iso_p7v6source.txt table. 


2 . 1 . Spatial Analysis 

To study the morphological characteristics of the 7 -ray 
emission in the field of G5.7-0.1, events in the energy range 
from 2 to 200 GeV were used. Also, only front-selected pho¬ 
tons, those converted in the first 12 scintillator planes of the 
Fermi-LAI Tracker, were included in the analysis to ensure 
better quality photon data and improved angular resolution. 
This gives the analyzed events a 68 % containment radius an¬ 
gle for normal incidence photons of < 0.3°. 

The gtbin tool was used to make counts maps of a 20° x 
20° region centered on the region of interest (ROI) defined 
during the event selection steps. Then the Galactic diffuse 
and isotropic extragalactic backgrounds were modeled and 
a 1° x 1° test statistic (TS) map was constructed using the 
gttsmap tool to determine the significance of 7 -ray emission 

5 The Science Tools package v9r31pl and support documents are dis¬ 
tributed by the Fermi Science Support Center and can be accessed at 
http://fermi.gsfc.nasa.gov/ssc 
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Figure 1 . The left panel shows the counts map of front converted events in the energy range from 2 to 200 GeV nearby SNR G5.7-0.1. Units are reported in 10 4 
counts deg - - with pixel binning 0.01°, and the map has been smoothed through convolution using a Gaussian distribution of width 0.2°. The left panel is a 1° 
X 1° image centered on G5.7-0.1 where the magenta crosses indicate the positions of 2FGL cataloged sources in the ROI {Nolan et al.|20l 2} and the magenta 
ellipses indicate the 95% confidence limits on the position of both. The right p anel shows the corresponding TS map created for the region. The green contours 
show radio emission observed using the VLA at 20cm (1.4 GHz) {Brogan et al.|2006{. The black cross indicates the position of detected OH maser emission 
l Frail et al. 1996 Hewitt et al. 2009) an d the magenta cross within the green contours shows the position of the coincident 7 -ray source 2FGL J1758.8-2402c 
and the V HE source 1 11 . S S J1800-240C { Aharon ian et al. 20081 Nolan et al. 2012 1 . The white curves in both panels show test statistic contours made using the 
TS map presented in the right panel. They correspond to detection significance measurements of 8 . 9, and 10cr. 


detection and precisely establish the position, and possible ex¬ 
tent, of a corresponding source. The gttsmap tool measures 
the test statistic across a defined grid of pixel positions. The 
values measured by the tool correspond to the logarithmic ra¬ 
tio of the likelihood of a point source existing in addition to 
the background model at a given point in the grid and the like¬ 
lihood of the model existing without an additional source at 
that grid position, 21og (L s+ b/Lb). 

The ROI for SNR G5.7-0.1 was centered at (072000 = 
269.73, Sj 2 ooo = -24.06) for the analysis. The left panel of 
Figure [I] shows the 1° x 1° countsmap made using the gtbin 
tool for the 2 - 200 GeV energy band, with a pixel binning 
of 0.01°. This map have also been smoothed using a Gaus¬ 
sian distribution of width 0.2°. Overlaid on the left panel are 
test statistic (TS) contours made using the TS map created 
with the gttsmap tool and shown in the right panel of the fig¬ 
ure. To create the TS map, four sources of 7 -ray emission 
were incorporated into the background source model. These 
included the two sources of diffuse 7 -ray background, both 
the Galactic and extragalactic components, the 7 -ray emis¬ 
sion from a nearby 7 -ray pulsar, PSR J1800-240, which is 
coincident with both 2FGL J1800.8-2400 and the VHE 7 - 
ray source HESS J1800-240B, and lastly the extended 7 -ray 
emission associated with W28, modeled in the 2FGL catalog 
by source J1801.3-2326e, which dominates the emission in 
the northeast corner of the countsmap. 2FGL J1800.8-2400 
was modeled as a point source with a power law spectral dis¬ 
tribution while 2FGL J1801.3-2326e was modeled as a disk 
with a log parabola distribution and a 0.39° radial extension, 
hist described by (|Abdo et al.|2010Jl. The parameters for both 
sources were fixed during the analysis at the values reported in 
the 24-month Fermi-LAT Second Source Catalog (Nolan et al. 

T 


2012 1 The vertical scale shows the measured TS values 


at 


each pixel position in the map, where \/TS is approximately 
the significance in units of the standard deviation. Using these 


6 The data for the 1873 sources in the Fermi LAT Second Source 
Catalog is made available by the Fermi Science Support Center at 
http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ 2 yr_catalog/ 


values, contours were made that represent significant emis¬ 
sion at 8 , 9, and 10er, given in white. 

Both panels also show a black cross representing the lo¬ 
cation of OH(1720 MHz) maser detection and overlaid radio 
contours, in green, from a 1.4 GHz VLA observation of G5.7- 
0.1 (Brogan et af||2006Jl. The location of 2FGL J 1758.8— 
2402c is coincident with the position of VHE 7 -ray source, 
HESS J1800-240C, and is marked by the magenta cross in¬ 
terior to the green radio contours (Aharoni an et al.| [2008). 
Finally, the magenta ellipses in both panels show the 95% 
confidence limit on the position of the two 2FGL 7 -ray point 
sources mentioned above (|Nolan et al.|2012]l. 


2 . 2 . Spectral Analysis 

To study the 7 -ray spectral energy distribution (SED) of 
G5.7-0.1, likelihood analysis was performed using both front 
and back converted events for the energy range 0.2-204.8 
GeV. The lower energy value is selected to avoid the highly 
variable effective area of the instrument at low energies and 
because of the large uncertainty of the Galactic diffuse model 
for energies below 0.2 GeV. The gtlike tool is used to model 
the 7 -ray flux for logarithmically spaced energy bins and to 
estimate the best-fit spectral parameters through likelihood 
maximization fitting. This requires modeling both the Galac¬ 
tic diffuse and isotropic extragalactic background and all 7 - 
ray point sources within 25° of the source of interest. The 
point sources are modeled in accordance with the 2FGL cat¬ 
alogue (Nolan et al.| |2012| ). The parameters of the galactic 
and extragalactic background models and all 2FGL sources 
within 5° of the source of interest were allowed to vary dur¬ 
ing the fitting. This included the parameters of the extended 
7 -ray source 2FGL J18 01.3- 2326e which was modeled in the 
same way as in Section [271] The modeled flux values result¬ 
ing from the likelihood analysis are reported in Figure [2] The 
statistically significant flux points represent the energy range 
0.2-25.6 GeV. The black error bars represent statistical uncer¬ 
tainties associated with the likelihood fitting method. The red 
error bars represent the systematic uncertainties. For the IRFs 
used in this analysis, the systematic uncertainties associated 
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with the effective area is energy dependent. For 100 MeV 
the uncertainty is 10%, decreasing to 5% at 560 MeV, and in¬ 
creasing to 20% at 10 GeV (Abd o et al.|2009a Nolan et al. 


|2012| and references therein). Additionally, the uncertainty 
in the flux of the Galactic diffuse background is also consid¬ 
ered. This was done by artificially varying the normalization 
of the Galactic background model from the best-fit flux val¬ 
ues by ± 6 % for each energy bin. Using the artificially ad¬ 
justed background values and freezing them, we remeasured 
the flux values and used the new measurements to approxi¬ 
mate the upper and lower limits of the uncertainty associated 
with the Galacti c diffuse backgroun d . This treatment is sim - 
ilarly applied inlAbdo et al.|(|2009c)>; |Castro & Slane|(|2010|); 
and |Castro et al. |2013[ >. 



Energy (MeV) 

Figure 2. Logarithmic plot of the high energy 7 -ray SED measured for SNR 
G5.7-0.1. The data points (black circles) represent flux values produced by 
the likelihood maximization process for energy ranges with >5cr detections, 
TS >25. The blue errorbars indicate statistical uncertainties resulting from 
the likelihood fitting while the red errorbars are systematic errors. The in¬ 
verted blue triangles correspond to calculated upper limits set on the energy 
flux for energy ranges with TS < 25. Included is a power law fit to the signif¬ 
icant data (red line) giving a measure of the spectral index, T = 2.27 zb 0.12. 


For the energy ranges with measured TS values < 25, the 
recorded 7 -ray flux has a significance of less than 5er. In 
these cases, the flux values are not reported but instead re¬ 
placed by an upper limit measured using the python based 
likelihood analysis tool, pyLikelihood. These flux upper lim¬ 
its are presented in Figure [2] alongside the flux measurements 
as inverted triangles. Also included is a y 2 -fit power law 
distribution for the seven significant flux values. The fitting 
measured a spectral index of T = 2.27 ±0.12. Integrating 
this flux model over the energy range, an estimate of the 
7 -ray flux for the energy range between 100 MeV and 100 
GeV measured at F 1 « 6.46 x 10 “' 1 erg/cnr/s. Using both 
kinematic distance measurements indicated by the adjacent 
OH(1720 MHz) maser emission velocity, 3.1 and 13.7 kpc, 
the luminosity over this band for each distance is estimated at 
L 7 31 ss 7.4 x 10 34 erg/s and L 7j 13.7 « 1.45 x 10 36 erg/s. 

2.3. XMM-Newton Analysis 

G5.7-0.1 was observed by XMM-Newton for ~ 32 ks 
on 16 September 2012 (ObsID 0691990101). The MOS 
and pn cameras were operated in full frame mode, and the 
"medium" optical filter was used. The data were reduced 
with the XMM-SAS software package (version 12.0.1) with 


calibration files current through XMM-CCF-REL-296, using 
the standard procedure for reducing XMM-Newton data out¬ 
lined in the XMM-Newton ABC Guide and the Birmingham 
XMM-Newton Guide. 

The MOS2 image of G5.7-0.1 is shown in Figure[3](r/g/!f). 
The image has been filtered to the energy range 0.5 ^2.6 keV 
and smoothed by a 36 arcsec Gaussian. Distinct arc-like arti¬ 
facts can be seen in the lower half of the image. These result 
from stray light associated with singly-reflected X-rays from 
the bright off-axis low mass X-ray binary GX 5-1. The arti¬ 
facts are exceptionally bright at higher energies, but modest 
below 2.6 keV due to the high column density for GX 5-1. 
Contours are from the 1.4 GHz VLA image shown in the left 
panel of Figure [3j( Brogan et al.|2006jl. 

We extracted spectra from the northern “stripe-free” half of 
G5.7-0.1 and found that the emission is adequately described 
by an absorbed solar-metallicity plasma with kT ~ 0.65 keV 
and Ah = (1.3 ±0.1) x 10 22 cm 2 , presumably associated with 
the SNR, accompanied by an absorbed bremsstrahlung com¬ 
ponent with kT = 7.2 keV and N H = 2.5 x 10 22 cm~ 2 , c onsis- 
tent with scattered emission from GX 5-1 ( jSmith et al.|20061 ). 
The contribution from the latter component prevents detection 
of any potential nonthermal emission from G5.7-0.1. 

As noted in Section [T] OH(1720 MHz) maser measure¬ 
ments place G5.7-0.1 at a distance of either ~ 3.1 kpc or 
~ 13 .7 kpc. The estimated distance to GX 5-1 is ~ 8 kpc 
((Smith et al.|2006|. The lower column density for G5.7-0.1 
indicates that it is closer than GX 5—1, suggesting that the 
distance ambiguity is broken in favor of the near distance. 

3. DISCUSSION 

Based in part on the observation of nonthermal X-ray and 
7 -ray emission, SNRs are known to accelerate particles to 
very high energies. These accelerated particles produce high 
energy emission through either leptonic processes, including 
nonthermal bremsstrahlung and inverse Compton scattering 
of ambient photons by accelerated electrons, or through the 
decay of neutral pions resultant from the collision of acceler¬ 
ated hadrons with ambient protons. To better understand the 
nature and origin of the broadband, nonthermal emission in 
the direction of G5.7-0.1, it was necessary to model it using 
methods described in depth by Castro et al. ( 2013) >. 


3.1. Modeling broadband emission 


Models simulating emission due to shock acceleration of 
electrons and protons were used to match the broadband SED 
of G5.7-0.1 particularly the observed 7 -ray emission. To an¬ 
alyze the 7 -ray SED of the accelerated particles, and assum¬ 
ing isotropy of both particle species, proton or electron, the 
spectral momentum distribution of the accelerated particles, 
dNi/dp, was modeled using a power law with an exponential 
cutoff. 


dN, 

, = exp 

dp 



( 1 ) 


where i denotes the particle species, and po.i is the correspond¬ 
ing exponential cutoff momentum for that species ( Castro &] 
Slan e|20f0j >. Transforming the particle distribution to energy 
space allows the exponential cutoff to be defined by an input 
energy, E 0l . The parameter a is the particle momentum dis¬ 
tribution index. The spectral index is related to a as F = a—2. 
For the case of re-acceleration of ambient cosmic rays, the 
index values for the electrons and protons differ at the outset, 
but reach the same values as the process approaches the steady 
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Figure 3. Radio (left) and X-ray (right) images of G5.7-0.1. Contours correspond to the 1.4 GHz VLA intensity (using a linear scale from 1 to 20 mJy beam -1 ). 
The arc-like stripes in the lower half of the X-ray image are stray-light artifacts produced by singly-reflected X-rays from the bright off-axis source GX 5—1, 
located to the southeast of the SNR (and outside of the telescope field of view). 


we assume a p = a e . The 


state (Tang & Chevalier|2015 i. Here 
ratio of the normalization constants of these particle distribu¬ 
tions, cij, establishes the electron to proton ratio, k ep = a e /a p , a 
value of particular interest to this analysis. These coefficients 
are implicitly determined by defining k ep during the modeling 
and setting the total integrated energy converted into particle 
acceleration within the SNR shell E C r = OEsn, where 6 rep¬ 
resents the efficiency of converting the total SN energy into 
shock acceleration of cosmic rays. The total mechanical en¬ 
ergy for the SN was set to the canonical value, E SN = 10 51 erg. 
The model assumes that the SNR is a spherical shell of shock- 
compressed material, with thickness R /12 and compression 


ratio 4, expanding into a uniform medium of density no (Cas¬ 
tro et al.|2013| >. This ratio is expected to be somewhat larger 


if efficient particle accelerating is in fa ct occurring ( | Warren | 
et al. 2005} Cassam-Chenai et al .{2008) >. 

Modeling of the 7 r u -decay emission mechanism is based on 
predictions fromjKamae et al.| <j2006]), while the synchrotron 
and inverse Compton scatte ring (IC) emission c omponents 
utilize models described by Baring et al. ( jl999| l. The in¬ 
verse Compton scattering model assumes an ambient pho¬ 
ton field of Cosmic Microwave Background (CMB) photons 
with kT CMB = 2.725 K. Finally, the nonthermal bremsstrahlung 
emission is modeled using the work of]Bykov et al.|(|2000(l. 
Further descrip tion of the m odel used in this work can be 
found in Castro et al.|( |2013j i and references therein. 

Several parameters were adjusted to model the observed 
broadband spectrum for both the leptonic and hadronic sce¬ 
narios. These included the momentum distribution index, 
a, with 4.00 representing a flat spectral shape, the post¬ 
shock magnetic field strength, ZB, the ambient density, no, 
and the cutoff energies, /'o,. An adjustment of Bi affects 
the flux of the radio synchrotron model. The nonthermal 
bremsstrahlung and 7r°-decay components of the high energy 
model are directly influenced by changing the value of hq. 
During the modeling, the cutoff energies were constrained to 
Eo,e < Eq-P to account for the differing radiative loss scenar¬ 
ios of accelerated electrons and protons. 

The electron to proton ratio has been measured by observa¬ 


tions of cosmic rays (CRs) at Earth’s surface to be k ep = 0.01 
( |Hillas||2005| >. There is also a known difficulty in injecting 
electrons into the DSA process. This means it is quite rea¬ 
sonable to expect k ep < 1 and potentiall y k ep << 1 at the par- 
ticle acceleration site ( |Park et al.|2015) . |Castro et al.| ( |2013) 
demonstrates that at k ep = 0.01 the hadronic component will 
dominate over the IC emission component for ambient den¬ 
sities larger than ~ 1 cm -3 over the energy range used to ob¬ 
serve G5.7-0.1 in this study. It was shown that the model 
predicts hadronic emission to dominate over the nonthermal 
bremsstrahlung for all of the particle distributions considered 
by Castro et al. ( 2013| as well. Therefore, the ratio was 
loosely constrained around k ep = 0.01 during the 7 r°-decay 
dominated modeling. It was also shown that an electron- 
proton ratio larger than k ep ~ 0.1 is required for the nonther¬ 
mal bremsstrahlung component to dominate the 7 -ray emis¬ 
sion in the Fermi -LAT band. Because of these considerations, 
nonthermal bremsstrahlung is not expected to dominate the 
applied models ([Castro et al.|2013jl. 

The final constraint fixed the total energy of the accelerated 
particles in the hadronic models to Eqr = Eqr, p +Ecr , e = 4 x 
10 50 erg. This assumes that 40% of the SNR energy from 
a 10 51 erg supernova is deposited into relativistic particles. 
This represents a reasonable upper limit on this energy and 
therefore sets a lower limit on the ambient density, n 0 , which, 
as mentioned previously, is constrained during the modeling. 
This should allow comparison between the estimated density 
predicted by the modeling with expectations for shocked ISM. 
The total CR energy and k ep are varied in the inverse Compton 
dominated scenario to fit the measured flux while limiting the 
contribution of the hadronic component over the Fermi-LAT 
energy range (Castro et al .|2013| ). 

3.2. SNR G5.7-0.1 emission 

Due to the ambiguity in the measured distance to G5.7-0.1, 
broadband modeling of the 7 -ray SED and radio observations 
of the SNR was conducted for both potential distances re¬ 
ported by Hewitt et al.]([2009[). Figure [4] shows the broadband 
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Table 1 

Inputted Model Parameters 


d 

Model 

a 

E o, p 

E(),e 

kep 

n 0 

b 2 

E C r,„ 

EcR,e 

(kpc) 



(TeV) 

(TeV) 


(cm -3 ) 

(jjG) 

(10 50 erg) 

(10 50 erg) 

13.7 

Hadronic 

4.00 

2.0 

2.0 

0.01 

22.0 

17.0 

3.894 

0.116 


Leptonic 

3.95 

15.0 

1.0 

0.5 

2.0 

1.5 

4.702 

2.298 

3.1 

Hadronic 

4.00 

2.4 

2.4 

0.001 

1.2 

12.0 

3.993 

0.007 


Leptonic 

3.95 

40.0 

1.2 

0.035 

0.22 

1.5 

3.881 

0.118 




Energy (MeV) 


Figur e 4. Broadband fits to radio observations (blue circles) fBrogan et al. 
[2006} . and Fermi -LAT observations (black circles) of SNR G.V/O.I with 
SNR distance assumed at 13.7 kpc (top) and 3.1 kpc (bottom). Flux points 
for the VHE 7-ra y sou rce HESS J1800-240C are included as green circles 
(Aharonian et al. [2008). The solid black curve represents the total modeled 
nonthermal emission. Also shown are the constituent spectra from IC (ma¬ 


genta; dot-dashed), 7r°-decay (blue; dashed), and nonthermal bremsstrahlung 
(green; dotted) emission. The lower energy component modeled the radio 
synchrotron emission. 


fits for both hadronically and leptonically dominated emission 
scenarios for the SNR being positioned at 13.7 kpc (top) and 
3.1 kpc (bottom) . The radio emission was measured by |Bro-| 
gan et al. (12006)1 using VLA observations of the SNR at 90 
cm (333 MHz), 20 cm (1.49 GHz) and 11 cm (272 MHz). 
These values and the Fermi-LAV flux points were fit by both 
the 7r°-decay and IC-dominated emission models described in 
the previous section. The flux of VHE 7 -ray source HESS 
J1800-240C, believed to be associated with G5.7-0.1, is in¬ 
cluded as green circles (Aharonian et al.|2008|. The results 


of the model fitting are shown in Table [I] Figure [5] focuses 
on the highest energy range of the spectrum in order to bet¬ 
ter understand the differences between the resulting fit of the 
two models at each distance. Also included in Figure[5]is the 
simple power law model fit applied previously and shown in 
Figure 2) 

Applying a model dominated by the hadronic emission 
mechanism to the 7 -ray emission in the direction of SNR 
G5.7-0.1 was favored at both distances. At a distance of 
13.7 kpc, the hadronic model used the electron-proton ra¬ 
tio observed on Earth, k ep = 0.01 with a spectral index of 
T = 2 . 0 , and a total energy of accelerated particles set to 
4 x 10 5l) erg. Additional parameters were adjusted until the 
model closely approximated the radio and 7 -ray flux mea¬ 
surements, including the HESS observations, and respected 
the calculated upper limits. Modeling indicated cutoff ener¬ 
gies equal to Eq.i = 2.0 TeV and an ambient density estimate 
of 770 , 13.7 = 22.0 cm -3 . Also a magnetic field of 17.0 /iG was 
required to accurately model the radio synchrotron flux. The 
resulting estimation of the energy converted into the acceler¬ 
ating cosmic ray protons (CRPs) was ~ 3.89 x 10 50 erg, with 
the remaining ~ 0.11 x 10 50 erg used in cosmic ray electron 
(CRE) acceleration. These results are consistent with previ¬ 
ous studies of SNR-MC systems that confirmed a hadronically 
dominated emission mechanism (Castro et al.|2013[|Auchettl 
let al.|2014| > 

The leptonic model at a distance of 13.7 kpc also closely 
approximated the measured data points and upper limits val¬ 
ues. However, the model required an electron-proton ra¬ 
tio of k ep = 0.5, significantly larger than that suggested by 
ground observation. Using this ratio and a spectral index of 
T = 1.95, different cutoff energies for both particle species 
were required to model the SED. Below ~ 20 GeV the lep¬ 
tonic components dominated the spectrum while above this 
energy range, the pion-decay channel was extended in order 
to better model the HESS VHE detection. The cutoff ener¬ 
gies for the electron and proton components were E () e = 1.0 
TeV and Eq p = 15.0 TeV respectively. An ambient density of 
no, 13.7 = 2.0 cm -3 , and a magnetic field of 1.5 /iG were also re¬ 
quired to adequately approximate the SED. The data required 
an increase in the total supernova energy from the canonical 
value resulting in an estimated energy used to accelerate CRPs 
of ~ 4.7 x 10 50 erg and ~ 1.5 x 10 50 erg to accelerate CREs. 
These are somewhat higher values that would require an SN 
having both efficient particle acceleration and slightly more 
energy than is typically expected. Several factors including a 
relatively high proton cutoff energy, a slightly higher required 
Esn, and a notably high required electron-proton ratio, pro¬ 
vide little support favoring a leptonically dominated emission 
scenario. 

Applying the hadronically dominated emission model to 
the broadband SED while setting an SNR distance of 3.1 
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Energy (MeV) Energy (MeV) 



Energy (MeV) Energy (MeV) 


Figure 5. Results of both hadronic (left) and leptonic (right) model fitting of Fermi-LXT detected 7-ray emission (black circles) and calculated upper limits 
(blue triang les) w ith SNR distance at 13.7 kpc (top) and 3.1 kpc (bottom). Flux points for the VHE 7-ray source HESS J1800-240C are also included as green 
circles lAharonian et al.|2008f. The red errorbars are systematic uncertainties and the black errorbars are the statistical uncertainties. The red line represents a 
power law model fit to the data points covering the range 0.2-25.6 GeV. The solid black represents the total nonthermal emission as predicted by the model. The 
constituent modeled spectra from inverse Compton emission (magenta; dot-dashed), 7r°-decay (blue; dashed), and nonthermal bremsstrahlung (green; dotted), 
are also shown. 


kpc required an electron-proton ratio of k ep = 0 . 001 , an order 
of magnitude lower than that observed on Earth, and a hard 
spectral index of T = 2.0. A lower k,, p such as this has also 
been seen in the modeling of the 7 -ray SEDs of SNRs like 
RXJ1713.7-3746 and W41 ( [Ellison et akpOlOKCastFo et al.| 
2013) >. This value resulted in a total energy deposited in CRP 
acceleration of ~ 3.97 x 10 50 erg and ~ 0.03 x 10 50 erg used 
to accelerate CREs. A magnetic field of IE = 4.0 /iG was re¬ 
quired to approximate the radio synchrotron data. An impor¬ 
tant difference between the model parameters for the hadronic 
scenarios at the two distances were the required ambient den¬ 
sities. For an SNR distance of 3.1 kpc the model required 
a value of no, 3.1 = 10 cm -3 , an order of magnitude smaller 
than that required for the 13.7 kpc model. These parameters 
are all quite reasonable and provide convincing evidence of a 
hadronic origin to the observed 7 -ray emission. 

The leptonically dominated model also fit the broadband 
SED adequately. The model used a spectral index of T = 1.95 
and required an electron to proton ratio of k ep = 0.035. This 
somewhat higher k ep was required to keep the IC contribu¬ 
tion the dominant component of the total 7 -ray flux below 
~ 20 GeV. Again two separate cutoff energies were required 
to accurately model the HESS VHE detection. These cut¬ 


off energies were To.,, = 1.2 TeV and Eq p = 40.0 TeV re¬ 
spectively. Setting an ambient density of «o, 3 .i = 0.22 and a 
magnetic field of 1.5 //G resulted in a strong fit to the SED. 
These parameters resulted in an estimated CRP acceleration 
energy of ~ 3.88 x 10 50 erg and ~ 0.12 x 10 50 erg converted 
to CRE acceleration. It is important to point out that while the 
leptonic component dominates the model, the required pion- 
decay component is quite significant. Therefore, it appears 
that at this SNR distance a significant leptonic component to 
the 7 -ray emission cannot completely be ruled out. However 
in both the hadronic and leptonically dominated modeling a 
significant hadronic component was required. 

Further analysis was carried out to see if additional photon 
fields could bolster the IC component enough that a strong 
fit to the data was possible while using an electron-proton ra¬ 
tio of k ep = 0.01 and a total energy in accelerated particles of 
Ecr = 4 x 10 50 erg. A photon field of temperature T = 25K, 
which is characteristic of interstellar dust, was used in this 
modeling. For an SNR distance of 13.7 kpc, an additional 
photon field with a density ~ 40 times that of the CMB was 
required. This is a significantly higher density than what is 
expected from a similar Galactic photon field (Strong et al. 
2000|l. At 3.1 kpc an additional photon field ~ 5.1 times the 



































































8 


Joubert, Castro, Slane & Gelfand 


density of the CMB was necessary to fit the SED. This re¬ 
sulted in a strong fit with no, 3.1 = 0.43 cm -3 , Bi = 2.5 /rG and 
reasonable acceleration energies, Ecr. p = 3.94 x 10 5() erg and 
EcR.e = 0.06 x 10 5() erg. This photon density is again too large 
for the region, leaving the hadronic dominated model the most 
convincing fit to the data at either distance. 

While a hadronically dominated model appear to be the 
most appropriate fit to the broadband SED at either distance, 
the lower column density measured during the XMM-Newton 


analysis of G5.7-0.1, described in Section 2.3 makes a strong 
argument in favor of the 3.1 kpc distance. The parameter val¬ 
ues required to model the SED at 3.1 kpc match closely to 
previous studies of similar SNR-MC systems that used this 


model (Castro & Slane|2010; |Slane et al. 

2012 

|Castro et al. 

2012j|Castro et al. 2013l|Ackermann et al. 

2013 

. Therefore it 


and significant SNR-MC 7 -ray source, located at 3.1 kpc. 


4. SUMMARY 

The work presented herein studied the spatial and spec¬ 
tral characteristics of the 7 -ray emission from the region sur¬ 
rounding the SNR G5.7-0.1 that is known to be interacting 
with nearby molecular clouds. This study found a significant 
region of 7 -ray emission in the direction of the SNR at lOrr 
for the energy range 2-200 GeV with the Fermi-LAI. We also 
successfully measured the 7 -ray emission spectrum for the 
0.2-204.8 GeV band, detecting significant flux values for the 
energy range 0.2-25.6 GeV. The radio and 7 -ray broadband 
emission spectrum was then fit to a model for emission from 
a cut-off power law distribution of relativistic protons and 
electrons. This modeling allowed conclusions to be drawn 
about the origin of the observed emission. In an effort to set 
constraints on the position of the SNR, the model was ap¬ 
plied while assuming two different potential SNR distances. 
In both cases the broadband SED was better explained by a 
model that assumed a hadronic origin of the 7 -ray flux than 
one assuming a leptonic dominated scenario. XMM-Newton 
observations also measured a column density for the SNR that 
supported an SNR distance of ~ 3 kpc. This potentially estab¬ 
lishes SNR G5.7-0.1 as an SNR-MC system independent of 
the SNR W28 cloud complexes and a significant 7 -ray source. 
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